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ABSTRACT 

The  spatial  distribution  of  surface  clutter  in  general  depends  on  radar  resolution, 
grazing  angle  and  scatterers  on  the  surface.  The  Weibull,  K-  and  lognormal 
distributions  are  commonly  used  to  approximate  the  clutter  spatial  distribution. 
Comparisons  among  all  three  distributions  are  reviewed  and  extended.  Statistical 
properties  of  the  Weibull  and  K-  distributions  in  the  log  domain  are  derived  and  then 
used  in  new  approaches,  named  as  unbiased  estimation  schemes,  for  faster  parameter 
estimation  for  the  Weibull  and  K-  distributions.  The  maximum  likelihood  estimates  for 
the  K-distribution  are  also  derived  without  the  need  of  the  derivative  of  the  Bessel 
funtion.  The  proposed  unbiased  estimation  schemes  provide  nearly  as  identical 
estimates  as  the  maximum  likelihood  scheme  does,  according  to  real  aperture  radar 
data  and  synthetic  aperture  radar  data  analysed. 
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Clutter  Spatial  Distribution  and  New  Approaches  of 
Parameter  Estimation  for  Weibull  and 
K-  Distributions 


EXECUTIVE  SUMMARY 

This  report  presents  the  result  of  a  study  of  the  clutter  spatial  distribution  carried  out 
under  the  task  AIR02/232  (AEW&C  SUPPORT). 

Clutter  return  from  sea  and  land  severely  limits  the  performance  of  radar  surveillance 
systems.  The  motivation  for  determining  the  correct  clutter  spatial  distribution  is  that 
the  long  'tails'  of  the  distribution  are  crucial  for  setting  target  detection  thresholds,  in 
order  to  minimise  both  false  alarms  and  missed  detections. 

The  clutter  spatial  distribution  is  complex  and  varies  from  situation  to  situation.  It 
depends  on  not  only  the  characteristics  of  random  scatterers,  but  also  radar  resolution 
and  the  geometry  of  illumination  (grazing  angle),  as  well  as  radar  frequency  and 
polarisation.  In  the  case  where  the  radar  resolution  is  low  and  the  iHuminated  surface 
is  homogeneous,  the  clutter  distribution  is  or  is  close  to  a  Rayleigh.  But  in  the  case 
where  the  radar  resolution  is  high  to  very  high  and  the  grazing  angle  is  low,  the  clutter 
distribution  of  even  a  homogeneous  land  or  sea  surface  significantly  diverges  from  a 
Rayleigh  due  to  the  effects  of  shadowing  and  multipath  propagation.  The  Weibull,  K- 
and  lognormal  distributions  are  commonly  used  to  represent  clutter  spatial 
distributions  for  various  situations.  However,  some  issues  related  to  these 
distributions,  such  as  faster  parameter  estimation  schemes  deserve  a  research  effort. 

The  contribution  of  this  report  is  three-fold. 

Comparisons  among  these  three  distributions  are  reviewed  and  extended.  It  is  found 
that  the  difference  between  the  Weibull  and  K-  distributions  is  small  when  the  shape 
parameter  is  within  a  certain  range  (often  clutter  distributions  are  well  within  this 
range).  Previous  comparisons  appearing  in  literature  are  found  to  be  incomplete. 

In  general,  clutter  distributions  are  often  plotted  on  a  semi-logarithm  scale,  i.e.,  clutter 
is  expressed  in  decibel  (dB).  Statistics  of  distributions  calculated  in  the  log  domain  are 
desirable.  The  mean,  standard  deviation  (or  variance)  and  so  forth  for  the  Weilbull  and 
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K-  distributions  in  the  log  domain,  which  were  previously  not  available,  are  derived.  In 
this  report  these  properties  are  later  used  in  the  parameter  estimation  schemes. 

Parameter  estimation  is  an  important  issue  and  the  maximum  likelihood  (ML) 
estimates  are  the  optimal.  Previous  research  concentrated  on  using  the  Lagrange 
method  to  find  the  ML  estimates.  This  seems  inappropriate  for  the  K-distribution  since 
the  derivative  of  the  embedded  Bessel  function  cannot  be  expressed  in  a  closed  form.  It 
is  therefore  very  difficult  to  find  the  ML  estimates  for  the  K-distribution  using  the 
Lagrange  method,  if  not  impossible.  This  report  proposes  a  direct  method  of  finding 
the  ML  estimates  for  the  K-distribution.  The  direct  method  does  not  require  the 
derivative  of  the  Bessel  function.  The  ML  schemes  often  serves  as  a  benchmark  for 
assessing  other  schemes. 

The  ML  method  provides  the  optimal  parameter  estimates.  Unfortunately,  except  for 
the  lognormal  distribution,  the  ML  estimates  for  both  the  Weibull  and  K-  distributions 
involve  iterative  computation  manipulating  sample  data,  which  may  not  be  practical  in 
real-time  applications.  In  real-time  data  processing,  fast  algorithms  are  desirable.  This 
report  proposes  two  unbiased  estimation  schemes  for  estimating  parameters  of  the 
Weibull  and  K-  distributions.  The  results  are  nearly  identical  to  the  ML  estimates  for 
various  data  sets  analysed,  but  the  computation  should  be  significantly  faster  than  that 
of  the  ML  algorithms. 

According  to  both  the  real  aperture  radar  data  and  synthetic  aperture  radar  data 
analysed,  the  K-distribution  seems  marginally  better  than  the  Weibull  distribution  in 
describing  the  land  clutter  distribution,  especially  at  the  high  end  of  the  distribution. 
Among  the  three  distributions,  the  lognormal  distribution  is  found  to  be  the  best  fit  for 
the  distribution  of  low  grazing  angle  (close  to  zero  degrees)  and  high  resolution  (the 
range  resolution  less  than  lm)  sea  clutter  data. 
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1.  Introduction 


Clutter  return  from  sea  and  land  severely  limits  the  performance  of  radar  surveillance 
systems.  The  motivation  for  determining  the  correct  clutter  spatial  distribution  is  that 
the  long  'tails'  of  the  distribution  are  crucial  for  setting  target  detection  threshold,  in 
order  to  minimise  both  false  alarms  and  missed  detections. 

Clutter  statistics  are  complex,  and  its  variations  vary  from  situation  to  situation.  It  is 
well  known  that  the  fluctuations  of  thermal  noise,  or  clutter  return  from  scatterers  of 
approximately  equal  magnitude  randomly  distributed  in  space,  obey  Rayleigh 
statistics.  However  it  has  been  found  that  when  measuring  a  complex  target  or  clutter, 
the  probability  of  relatively  large  signal  occurrence  is  higher  than  that  predicted  from  a 
Rayleigh  model.  Such  distributions  are  said  to  have  long  tails.  In  general,  if  radar  looks 
at  natural  land  with  a  low  resolution,  the  clutter  distribution  is  or  is  very  close  to 
Rayleigh.  However  in  the  case  of  high  and  very  high  resolution  radars,  the  resultant 
distribution  is  generally  much  wider.  Additionally,  clutter  return  at  low  grazing  angles 
does  not  obey  the  Rayleigh  distribution  because  of  effects  of  shadowing  and  multipath 
propagation.  Therefore,  the  Weibull,  K-  and  lognormal  distributions  are  most 
commonly  used  to  describe  clutter  statistics  (the  Rayleigh  distribution  is  a  special  case 
of  the  Weibull  or  K-  distributions).  Characteristics  of  these  three  distributions  have 
been  considerably  studied  and  reported  (Jao,  1984,  Watts,  1985,  1987,  Oliver,  1993, 
Blacknell,  1994,  Antipov,  1998,  Long,  1992, 2001,  Billingsley,  2002). 

The  contribution  of  this  report  is  three-fold. 

Comparisons  among  these  three  distributions  are  reviewed  and  extended.  It  is  found 
that  the  difference  between  the  Weibull  and  K-  distributions  is  small  when  the  shape 
parameter  is  within  a  certain  range  (often  clutter  distributions  are  well  within  this 
range).  Previous  comparisons  such  as  those  given  by  Billingsley  (2002)  are  found  to  be 
incomplete. 

In  general,  clutter  distributions  are  often  plotted  on  a  semi-logarithm  scale,  i.e.,  clutter 
is  expressed  in  decibel  (dB).  Statistics  of  these  distributions  calculated  in  the  log 
domain  are  desirable.  The  mean,  standard  deviation  (variance)  and  so  forth  for  the 
Weilbull  and  K-  distributions  in  the  log  domain,  which  were  previously  not  available, 
are  derived.  In  this  report  these  properties  are  later  used  in  the  parameter  estimation 
schemes. 

The  main  contribution  of  the  report  is  relevant  to  parameter  estimation  for  both  the 
Weibull  and  K-  distributions.  Previous  research  concentrated  on  using  the  Lagrange 
method  to  find  the  maximum  likelihood  (ML)  estimates.  This  seems  inappropriate  for 
the  K-distribution  since  the  derivative  of  the  (imbedded)  Bessel  function  cannot  be 
expressed  in  a  closed  form.  It  is  therefore  very  difficult  to  find  the  ML  estimates  for  the 
K-distribution  using  the  Lagrange  method,  if  not  impossible.  This  report  proposes  a 


1 


DSTO-RR-0274 


direct  method  to  find  the  ML  estimates  for  the  K-distribution.  The  direct  method  does 
not  require  the  derivative  of  the  Bessel  function. 

Theoretically  the  ML  estimates  are  the  optimal.  Unfortunately,  iterative  algorithms 
manipulating  sample  data  are  involved  in  determining  the  ML  estimates  for  both  the 
Weibull  and  K-  distributions,  which  may  not  be  practical  for  real-time  applications.  In 
real-time  data  processing  fast  algorithms  are  desirable.  This  report  also  proposes  two 
unbiased  parameter  estimation  schemes  for  the  Weibull  and  K-  distributions.  The 
estimates  are  simple  and  yet  found  to  be  nearly  identical  to  the  ML  estimates  for  both 
the  real  aperture  radar  (RAR)  and  synthetic  aperture  radar  (SAR)  data  analysed. 

2.  Weibull,  K-  and  Lognormal  Distributions 

2.1  Weibull  Distribution 

The  Weibull  distribution  has  been  found  to  provide  a  good  fit  to  both  sea  and  land 
clutter  (Billingsley,  2002).  The  probability  density  function  (pdf)  of  the  Weibull 
distribution  on  a  linear  scale  is, 

p  (x)  =  bcxb~l  exp(-cxA )  x  >  0  (1) 

where  a -\lb  and  c  are  often  referred  to  as  the  Weibull  shape  and  scale  parameters, 
respectively. 

The  cumulative  distribution  function  (cdf)  of  the  Weibull  distribution  is, 

P(x)  =  1  -  exp(-cxA )  (2) 

Since  clutter  is  in  general  expressed  in  decibel  (dB),  often  we  may  prefer  statistics  of 
clutter  distribution  to  be  expressed  on  the  dB  scale.  In  this  report,  we  use  x  to  denote 
the  random  variable  on  the  linear  scale,  y  =  lnx  on  the  natural  logarithm  (NL)  scale  and 
z  =  io  log*  on  the  dB  scale.  Because  the  relationship  between  y  and  z  is  only  a  constant 

scale  factor,  i.e.,  2  =  ,  any  properties  of  either  one  would  also  apply  to  the  other. 

Into 

This  report  uses  both  y  and  z:  y  for  neat  mathematical  manipulations  and  z  to 
express  clutter  in  dB. 

Knowing  the  pdf  of  a  variable  x  and  a  unique  function  y  =  y(x) ,  we  want  to  find  the 
pdf  of  y  .  Because  the  interval  (x,x  +  dx)  is  mapped  onto  the  interval  (y,y  +  dy),  the 
probability  of  the  variables  falling  in  the  intervals  is  the  same.  That  is, 

p(x)dx=p(y)dy  (3) 
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Because  p(y )  has  to  be  >  0 ,  so  we  have 


p(y)  =  p(x)\ 


dx 

dy 


(4) 


In  the  NL  and  dB  domains,  respectively,  the  pdfs  of  the  Weibull  distribution  become, 
P(y)  =  bc{ey  f  exp|-  c(ey  J  j  (5) 

p(z)=j~b4 102/10  F  exp^-c^O2'10)^  (6) 


where  k0  =  10/ln  10. 

Table  1  lists  statistics  of  the  Weibull  distribution.  Properties  on  the  linear  scale  are 
available  in  literature  (e.g.,  Billingsley,  2002,  Long  2001).  Derivations  of  the  properties 
in  the  NL  and  dB  domains  are  given  in  Appendix  A. 


Table  1:  Statistics  of  the  Weibull  distribution  in  the  linear,  NL  and  dB  domains,  respectively. 


Linear  domain 

NL  domain 

DB  domain 

Pdf 

p(x)  =  bcxb~l 
exp  {-cxb  ) 

p{y)  =  bc[eySf 
exp^-c(ey)6  j 

p(z)=-Lbc[ 10z/1°)£' 
k0 

exp^-c^O2710)6  j 

Cdf 

P(x)  =  l-exp(-cx*) 

POO  =  l-expf-c(ey)f  j 

jj 

Mean 

r(l  +  a) 

c" 

— io'  +  lnc) 

_y(r  +  inc) 

Median 

— [ln(ln2)  -  Inc] 
b 

—  [ln(ln2)-lnc] 
b 

Mode 

U)  f°'a<l 

Inc 

b~ 

-k  — 

Ko  , 
b 

Variance 

ms 

x2 

7T2 

6b2 

l .  2  2 

k0  Tt 

6  b2 

Std 

fra+2fl)  ^ 

lr2(l  +  a)  J 

1/2 

X 

1,  -O 

^  r> 

ko7T_ 

■JZb 

where  y  =  -yr(0)  (1)  =  0.5772156649 ,  is  Euler's  Gamma  constant. 
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2.2  K-Distribution 

Many  researchers  have  reported  that  both  sea  clutter  and  land  clutter  obey  the  K- 
distribution  (Jao,  1984,  Watts,  1985  and  1987,  Oliver,  1993).  The  K-  distribution  has  been 
well  interpreted  by  Ward  et  al,  who  named  it  the  compound  K-distribution  (Ward  et  al, 
1990).  The  particular  significance  of  the  K-distribution  is  its  compound  nature,  which 
interprets  clutter  comprising  two  components:  a  fast  component  (correlation  time 
about  10ms  for  sea  clutter  at  X-band),  described  as  a  Rayleigh  distribution,  models  the 
local  speckle;  and  a  slow  component  (correlation  time  in  the  order  of  seconds  for  sea 
clutter  at  X-band),  characterised  by  a  gamma  distribution,  models  the  underlying  mean 
of  the  speckle  level,  also  referred  to  as  texture  (Gird  et  al,  2001,  Gird  and  Greco,  2002).  It 
has  been  found  that  the  compound  K-distribution  model  accurately  describes  real  sea 
clutter  collected  by  a  vertically  polarised  X-band  radar  at  a  low  grazing  angle  (Farina, 
2002). 

There  are  at  least  two  versions  of  the  K-  distribution,  one  for  square  law  detection 
(power  domain)  and  the  other  for  amplitude  detection  (voltage  domain).  However 
they  are  slightly  different.  The  expression  for  square  law  detection  given  by  Jao  (1984) 
is. 


a-1 


P(*)  = 


2  a 


r(ar)x 


ax\  2 


Ka-x 


_  ax  ,  A 
2,1 -  x  >  0 


(7) 


where  Ku(-)  is  the  modified  Bessel  function  of  the  second  kind. 

The  expression  for  amplitude  detection  is  (Watts  1985, 1987,  Antipov  1998,  Long  2001), 


P(t)  = 


2c 


T(a) 


2 


K-a-M 


(8) 


By  substituting  u  =  x1'2  and  d  =  ijaix  into  (7),  we  have. 


(9) 


On  the  surface,  (9)  appears  exactly  the  same  as  (8),  but  they  are  actually  not.  In  (8) 
parameters  a  and  c  are  assumed  to  be  independent  while  parameter  d  in  (9)  is 
dependent  on  parameter  a .  This  difference  will  lead  to  different  expressions  when  we 
manipulate  the  pdf  with  respect  to  parameters  (parameters  are  treated  as  variables), 

such  as  —p(x),  in  order  to  find  the  parameter  estimates  using  the  ML  method,  for 
8a 
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example.  However  in  terms  of  the  statistics  of  the  distribution  (parameters  are  treated 
as  constants),  both  (7)  and  (8)  should  be  the  same.  In  this  report  we  adopt  (7)  as  the  K- 
distribution,  which  has  been  widely  used  (Billingsley,  2002,  Oliver  1993),  consistent 
with  the  Weibull  and  lognormal  distributions  for  the  square  law  detection. 

The  pdfs  of  the  K-distribution  when  the  variable  is  expressed  in  the  NL  and  dB 
domains,  respectively,  are. 


(11) 


The  statistics  of  the  K-distribution  in  the  linear,  NL  and  dB  domains,  respectively,  are 
given  in  Table  2.  The  corresponding  derivations  can  be  found  in  Appendix  A. 


Table  2:  Statistics  of  K-distribution  in  the  linear ,  NL  and  dB  domains,  respectively. 


Linear  domain 


2 a  f ax i  2 
r(or)*l  x  , 


NL  domain 


0+1 

p(y 0  =  — f-eyl  2 


■Ka. i 


P  00  =  1- 


2  I  ae 


T(a)  {  x 


-y  +  ip(  ' (or)  —  In — 
x 


y/m{a)  +  ^~ 


yr()  (a)  +  - 


DB  domain 


piz)  =  -L-2-(£L  io^10 

k0  r(a)U 


■Ka- 1  2j"lOz/1° 


P(z)  =  1- 


2  f a\0y'w 


T(a)  7 
.  Ial0yn°) 


k0\  i//w(a)  +  - 


Where  v<n)(x)  is  the  Poly  gamma  function,  defined  as  the  nth  derivative  of  the  logarithm 


of  the  Gamma  function, - -  (in  T(;t)). 
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It  can  be  proven  that  a  Weibull  distribution  with  b  =  1  is  equivalent  to  a  K-distribution 
with  a  =  oo ,  and  both  distributions  are  identical  to  the  negative  exponential 
distribution,  i.e.,  the  Rayleigh  distribution,  as  (Jao,  1984), 


p(*)  =  -rexp(-x/x)  (12) 

x 

A  Weibull  distribution  with  b  =  1/2  is  again  identical  to  a  K-distribution  with  or  =  1/2 
(Davidson,  et  al,  2002).  This  can  be  easily  proven  if  we  note  (Gradshteyn  and  Ryzhik, 
1980), 

K1/2(x)  =  ^exp(-x)  (13) 

2.3  Lognormal  distribution 

The  distribution  of  a  random  variable  x  is  lognormal,  if  tax  is  normally  distributed 
with  mean  p  and  variance  s1 .  The  pdf  of  the  lognormal  distribution  is  given  by. 


p(x)  = 


1  1 


( 


•Jins  x 


—exp 


V  2s 


1  -(In  x-ju)2 


x  >0 


The  cdf  of  the  lognormal  distribution  is 


P(x)  =  ^[l  +  sign(u)erf{\  u  |)] 


(14) 


(15) 


where  u  =  ln*^;  signQ)  is  the  sign  function;  and  er/Q  is  the  error  function,  defined  as 
V2 s 


0  u 

erf  (m)=  —t=  J  exp(-X  2)dt  u>0 

Vxr  o 

(16) 

The  pdfs  of  the  lognormal  distribution  when  the  variable 
domains,  respectively,  are, 

is  expressed  in  the  NL  and  dB 

pW=V^e4l7(,,''')) 

(17) 

PW  t„V2?SeXP(  2A0  ) 

(18) 
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Compared  to  the  Weibull  and  K-  distributions,  the  convergence  of  the  cdf  of  the 
lognormal  distribution  is  the  slowest,  i.e.,  its  'tail'  is  longest.  Therefore,  the  lognormal 
distribution  has  been  used  extensively  in  cases  where  echo/clutter  fluctuations  attain 
large  values  for  much  higher  percentages  than  Weibull  and  K-distribution  statistics 
(Long,  2001). 

Table  3  lists  statistics  of  the  lognormal  distribution  in  the  linear,  NL  and  dB  domains, 
respectively. 


Table  3:  Statistics  of  the  lognormal  distribution  in  the  linear,  NL  and  dB  domains,  respectively. 


B'SBffBTSi 

Linear  domain 

NL  domain 

DB  domain 

Pdf 

,  ,  11 
pW=  t — 

fins  * 

(  1  2*) 

exp - r-(lnx  -//) 

2sl  ) 

PW~f2 xs 

P(z)  .  r~ — 

«0  V2 7tS 

(  1  \ 

1  ,  Z  .2 

exp 

\  is  *o  y 

Cdf 

P(x)  =  ^  [l  +  sign(u)erf  (|  u  |)] 

.  In  x-u 

where  u  =  — == — 

V2s 

p(y)  =  ^[l  +  sign(u)erfQ  u  |)] 

where  u  = 

V2s 

P(z )  =  \ll  +  sign(u)erfQ  u  |)] 

,  zlk0-p 

where  u  -  — ■=■ — - 
V2  s 

Mean 

exp (m+s2  / 2) 

k0y 

Median 

expCa) 

M 

k0ju 

Mode 

expC u-s2) 

M 

hP 

Variance 

[exp(52)-l]jc2 

s2 

h2s2 

Std 

[exp(s2)-l]I/25c 

s 

k0s 

3.  Comparisons  Among  Weibull,  K-  and  Lognormal 

Distributions 


Because  the  logarithmic  transform  is  non-linear  with  low  values  being  stretched  and 
high  values  compressed,  the  same  distribution  on  the  linear  scale  and  the  dB  scale 
might  look  significantly  different.  As  clutter  is  generally  expressed  in  dB,  we  mainly 
concentrate  our  discussion  of  distributions  on  the  dB  scale. 

3.1  Same  Mean  and  Same  Standard  Deviation  in  the  Linear  Domain 

We  assume  x  =  1,  as  it  is  only  a  matter  of  scaling  in  the  linear  domain  or  a  matter  of 
shift  in  the  dB  domain  if  the  actual  mean  value  is  different.  We  also  assume  std(x) 
(standard  deviation  of  x )  is  the  same  for  all  three  distributions  in  each  case  for  the 
purpose  of  comparison.  In  particular.  Table  4  shows  mean  and  std  values  when 
calculated  in  the  dB  domain  for  four  cases,  corresponding  to  a  =  1,1.5, 2, 3  for  the 
Weibull  distribution,  respectively.  A  Weibull  distribution  with  a  =  1  is  the  same  as  a  K- 
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distribution  with  «  =  <*>,  and  both  simplify  to  the  Rayleigh  distribution.  A  Weibull 
distribution  with  a  =  2  is  also  the  same  as  a  K-distribution  with  a  =  0.5  .  Although  in 
each  case,  E(x)  and  std(x)  are  the  same  for  all  three  distributions,  E(z)  and  std(z)  for 
the  three  distributions  diverse  remarkably  especially  when  std(x)  becomes  larger. 


Table  4:  Distribution  comparison  from  the  perspective  on  the  mean  and  std  in  the  dB  domain 
giving  the  mean  and  std  in  the  linear  domain  to  be  the  same  for  all  three  distributions. 


Weibull 

K- 

Lognormal 

Case  1:  E(x)  =  1;  a  =  1  (srd(x)  =  l) 

Parameters 

a  =  1;  c  =  1 

a  =  °° 

H  =  -0.3466;  5  =  0.8326 

E(z) 

-2.5dB 

-2.5dB 

-1.5dB 

std(z) 

5.6dB 

5.6dB 

3.6dB 

Case  2:  E(x)  =  1;  a  =  1.5  (std(x)  =  1.5477) 

Parameters 

a  =  1.5;  c  =  1.209 

a  =  1.4334 

=  -0.6112;  5  =  1.1056 

E(z) 

-5.0dB 

-4.2dB 

-2.7dB 

std  (z) 

8.4dB 

7.1dB 

4.8dB 

Case  3:  E(x)  =  1;  a  =  2  (std(x)  =  2.2361) 

Parameters 

a  -  2;  c  =  1.414 

o 

il 

£3 

p  =  -0.8959;  5  =  1.3386 

E(z) 

-8.0dB 

-8.0dB 

-3.9dB 

std(z) 

ll.ldB 

ll.ldB 

5.8dB 

Case  4:  E(x)  =  1;  a  =  3  ( st d(x)  =  2.2361) 

Parameters 

a  =  3;  c  =  1.8171 

«  =  0.1111 

p  =  —1.4979;  5  =  1.7308 

E{z) 

-15.3dB 

-33.8dB 

-6.5dB 

std  (z) 

16.7dB 

39.8dB 

7.5dB 

When  1  <  a  <  2 ,  the  Weibull  and  the  K-  distributions  are  different,  but  the  difference  is 
small.  This  will  become  clear  when  we  look  at  their  pdfs  and  cdfs  shortly.  In  this  range, 
std(z)  of  the  K-distribution  is  less  than  the  value  of  the  Weibull  distribution,  i.e.,  the 
shape  of  the  K-distribution  is  slightly  narrower  than  the  shape  of  the  Weibull 
distribution.  However  in  the  range  of  a >2,  the  shape  of  the  K-distribution  becomes 
broader  than  the  shape  of  the  Weibull  distribution.  The  shape  of  the  lognormal 
distribution  is  always  the  narrowest. 

Figure  1  shows  pdfs  of  Case  1,  2,  and  3,  respectively.  In  order  to  see  the  details  of  the 
distributions  at  the  low  probability  ends,  pdfs  on  the  log  scale  are  also  show  in  the 
figures.  It  is  easy  to  shown  that  the  pdf  of  the  lognormal  distribution,  when  expressed 
on  the  log  scale,  becomes  a  quadratic  function,  the  shape  of  the  pdf  is  therefore  a  face 
down  parabola  symmetrical  to  the  vertical  line  of  z  =  E(z) .  The  pdf  shapes  of  both 
Weibull  and  K-  distributions  are  also  similar  to  parabolas  (but  are  not  parabolas),  and 
their  quasi-symmetric  axes  are  not  vertical  but  tilted  some  degrees,  leading  to  smaller 
pdf  values  at  the  high  end.  The  Tail'  of  the  Weibull  and  the  k-  distributions  are 
therefore  not  as  long  as  the  'tail'  of  the  lognormal  distribution. 
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Figure  1:  Comparisons  of  pdfs  ofWeibull,  K-  and  lognormal:  (a)  a  =  1 ,  (b)  a  =  1.5  and  (c)  a  =  2 
respectively.  Pdfs  on  the  log  scale  show  details  of  distribution  in  the  low  probability  ends. 
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To  compare  the  convergence  of  cdf,  we  may  plot  cdfs  onto  the  Weibull  scale 

(Billingsley,  2002).  The  Weibull  scale  maps  the  cdf  P{z)  into  1 0  log^ln  ^—j,  so  the  cdf 

of  the  Weibull  distribution  becomes  a  straight  line  with  a  slope  of  b .  Figure  2  shows 
cdf  comparisons  among  the  three  distributions.  It  was  reported  that  among  the 
distributions,  the  K-  distribution  converges  fastest,  (curving  up),  the  Weibull  in  the 
middle  (straight  line)  and  the  lognormal,  the  slowest  (curving  down)  (Billingsley, 
2002).  From  Figure  2  we  can  see  that  this  conclusion  is  incomplete.  In  fact  only  when 
a  >2,  is  the  conclusion  correct,  while  in  the  range  of  l<a<2  the  cdf  of  the  K- 
distribution  also  curves  down  on  the  Weibull  scale.  Therefore  in  the  1  <  a  <  2  range,  the 
fastest  convergence  is  the  Weibull,  followed  by  the  K-  and  then  the  lognormal.  As 
mentioned  before,  when  a  =  1  or  2,  the  Weibull  and  the  K-  distributions  are  identical. 


Figure  2:  cdf  comparisons  using  the  Weibull  scale.  The  cases  a -1.5  and  a -2.5  are  shown,  in 
which  the  cdf  of  the  K-distribution  curves  down  and  up,  respectively,  indicating  that  the 
convergence  of  the  K-distribution  can  be  slower  or  faster  than  that  of  the  Weibull  distribution 
depending  on  the  corresponding  shape  parameter. 

Another  way  to  compare  the  cdfs  is  to  plot  log(l-P(z)),  so  the  details  of  the 
convergence  of  the  cdf  can  be  viewed  clearly  (It  is  understood  that  if  we  swap  the  axes, 
then  the  plot  will  become  the  so  called  quantile  distribution).  Figure  3  compares  the 
convergence  of  the  three  cdfs  with  a  =  1.5  and  a  =  2.5,  respectively.  The  lognormal 
distribution  converges  slowest,  and  requires  a  few  dB  more  than  the  others  in  order  to 
achieve  the  same  cdf  value  of  0.9999. 

In  practice,  the  clutter  distribution  in  most  cases  falls  in  the  range  1  <  a  <  2 ,  where  we 
have  seen  the  differences  between  the  Weibull  and  k-  distributions  are  very  small. 
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(a)  (b) 


Figure  3:  Cdf  comparison.  The  term  of  1  -P(z)  is  plotted  on  the  log  scale  to  shoxvn  the  details  of 
the  convergence  of  the  curves. 

3.2  Same  Median  and  Same  Standard  Deviation  in  dB  Domain 

It  can  be  seen  in  Table  4  that  for  the  same  std(x)  in  the  linear  domain,  the  lognormal 
distribution  usually  has  the  smallest  std{z)  in  the  dB  domain.  If  the  same  std  value  is 
given  in  the  dB  domain,  we  can  imagine  that  the  lognormal  distribution  has  even 
broader  shape.  The  convergence  of  its  cdf  is  then  even  slower,  i.e.,  its  'tail'  is  even 
longer. 

For  illustration  purposes.  Figure  4  shows  pdfs  in  which  we  assume  the  median  value  in 
the  dB  domain  to  be  OdB,  and  the  std  to  be  5.6dB  (corresponding  to  a  =  l  for  the 
Weibull  distribution)  and  ll.ldB  (corresponding  to  a  =  2  for  the  Weibull  distribution), 
respectively.  Compared  to  Figure  1,  where  the  same  mean  and  the  same  std  in  the 
linear  domain  is  assumed,  the  difference  between  the  lognormal  distributions  and  the 
other  two  widens.  As  shown  in  Figure  5,  for  the  lognormal  distribution  the  value  of  z 
has  to  be  more  than  lOdB  and  20dB  higher  than  that  of  the  Weibull  and  K-distributions, 
respectively,  to  achieve  the  same  cdf  value  of  0.9999.  This  demonstrates  the  importance 
of  the  clutter  distribution  studies. 


11 


DSTO-RR-0274 


(a) 


Figure  4:  Comparisons  of  pdfs  of  the  Weibull,  K-  and  lognormal  with  the  same  median  of  OdB, 
and  the  same  std  of  (a)  5.6dB  (equivalent  to  a  =  1  in  the  Weibull  distribution)  and  (b)  ll.ldB 
(equivalent  to  a  =  2  in  the  Weibull  distribution)  in  the  dB  domain. 


Figure  5:  Comparisons  ofcdfs  of  the  Weibull,  I 
and  the  same  std  of  (a)  5.6dB  (equivalent  to  a  ■■ 
(equivalent  to  a  =  2  in  the  Weibull  distribution) 


Z 


—  Weibull  &  K - Lognormal 


(b) 

and  lognormal  with  the  same  median  of  OdB 
=  1  in  the  Weibull  distribution)  and  (b)  ll.ldB 
in  the  dB  domain. 
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4.  Parameter  Estimation 


An  important  issue  in  the  clutter  distribution  studies  is  parameter  estimation.  For  a 
given  pdf,  its  parameters  determine  all  the  corresponding  statistics. 

4.1  Maximum  Likelihood  Method 

A  standard  approach  to  determining  parameters  is  the  ML  method,  which  provides 
optimal  parameter  estimates  in  the  sense  that  these  estimates  are  the  most  probable 
parameters  for  the  given  data  and  the  pdf,  if  there  is  no  prior  knowledge  (Oliver,  1993, 
Antipov,  1998). 

If  n  samples,  xl,x2,---,xn  are  given  to  estimate  m  parameters  0l,02,---,0m  of  a  given 
distribution  function,  then  the  joint  pdf  of  xl,x2,---,x„  is  the  product  of  the  marginal 
pdfs: 


f(x1,x2,---,x„\01,02,---,0m)  =  flf(xi;0i’02’---’0m) 

1=1 


(19) 


The  estimates  0l,02,---,0m  which  maximise  the  joint  pdf  (19)  are  referred  to  as  the  ML 
estimates.  In  general,  the  estimates  01,02,---,0m  can  be  found  using  the  Lagrange 
method. 


Since  the  logarithm  function  is  monotonic,  and  it  transfers  a  product  into  a  summation, 
the  ML  estimates  can  be  obtained  by. 


0O\  ;=i 

,^hnf{xrA,02-,0m)=o 

aU2  ,=i 


fin  f{xi-,0l,02,-,0m)  =  O 

i=i 


(20) 


4.1.1  ML  Estimates  of  Lognormal  Distribution 

Using  (20)  it  is  straightforward  to  derive  the  ML  estimates  for  the  lognormal 
distribution: 
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1  n 

ji  =  ~YMxl 

«/= l 

s 2  =-Z(lnjci  ->“)2 

«/= l 


(21) 


These  results  in  fact  are  the  mean  and  the  variance  of  the  samples  calculated  in  the  NL 
domain: 


J  M  =  {y)  =  E(y) 

|s2  =(>'2)-(y)2  =var(y) 


(22) 


where  (■)  denotes  the  ensemble  average. 

4.1.2  ML  Estimates  of  Weibull  Distribution 

Appling  (20)  to  the  Weibull  distribution,  it  is  not  difficult  to  derive  the  following  ML 
estimates  for  the  Weibull  distribution  (Oliver,  1993): 


c  n 


b 


.  (=i  i=i 


(23) 


Equation  (23)  in  general  cannot  be  solved  directly,  and  requires  an  iterative  solution. 
The  iteration  can  be  written  as. 


.  ,  .  V  *  =  1,2,-  (24) 

±-  =  ckE\xb*-'y  -£(v) 

h  K  J 

A 

Without  losing  generality,  we  can  let  b0  =  1 . 

4.1.3  ML  Estimation  for  K- Distribution 

Given  the  fact  that  the  derivative  with  respect  to  a  of  the  modified  Bessel  function  of 
the  second  kind  is  not  in  a  closed  form  (in  our  case,  a  is  not  only  the  part  of  the 
variable  but  also  the  order  of  the  function  making  the  derivative  even  more  complex), 
the  ML  estimates  of  the  K-distribution,  even  in  an  iterative  form,  are  difficult  to  find 
using  (20),  if  not  impossible.  Using  the  Mellin  transform,  Oliver  (1993)  has  derived  an 
asymptotic  expression  for  the  K-  distribution,  in  which  the  K-distribution  is  expressed 
as  a  Gamma  distribution  with  a  modification  factor.  The  corresponding  approximate 
ML  estimates  for  the  asymptotic  expression  have  then  been  derived.  However,  the 
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resultant  ML  estimates  is  only  applicable  to  the  limit  of  large  look  number  of  multi¬ 
look  data  (e.g.,  radar  measurements  are  averaged  in  the  power  domain,  such  as  the 
multi-look  SAR  data).  Therefore  the  formulae  do  not  suit  for  the  parameter  estimation 
for  (7)  in  which  the  look  number  is  assumed  to  be  one.  Joughin  et  al  (1993)  have 
proposed  a  direct  approach  to  find  the  ML  estimates  for  the  K-  distribution.  The 
method  expresses  the  K-distribution  as  a  combination  of  gamma  and  Rayleigh 
distributions  to  simplify  the  calculation,  but  the  estimates  have  to  be  numerically 
searched. 

In  fact,  if  the  calculation  is  not  required  in  real-time,  we  can  apply  an  iterative 
algorithm  searching  the  optimal  a  which  maximises  the  log-likelihood  function  (i.e., 
using  (19)  to  search  the  maximum  of  the  function  rather  than  using  (20)  to  search  the 
derivative  of  the  function  to  be  equal  to  zero),  because  the  direct  search  does  not 
require  the  derivates  of  the  Bessel  function.  Luckily,  such  a  direct  ML  estimates  search 
for  the  K-distribution  function  only  requires  a  one-dimensional  search.  The  parameter 
3c  does  not  need  the  optimal  search,  as  the  ML  estimation  of  3c  is  E(x) .  The  direct  ML 
estimates  search  scheme  for  the  K-  distribution  can  be  written  as  (omitting  constant 
terms). 


*=-£*/  =E(X) 
«/=l 


[2JO 

/  i= 1 

(25) 


4.2  Unbiased  Estimates 

Theoretically  the  ML  method  provides  the  optimal  parameters  which  are  the  most 
probable  parameters  for  a  given  distribution  and  sample  data.  In  the  asymptotic  limit, 
the  ML  estimates  produce  the  smallest  errors  that  may  be  calculated  via  the  Cramer- 
Rao  lower  bound  (Blacknell,  1994).  Unfortunately  except  for  the  lognormal 
distribution,  for  which  the  ML  estimates  can  be  directly  calculated  via  the  calculation 
of  the  mean  and  variance  (i.e.,  the  first  and  second  moment  of  the  data)  in  the  log 
domain,  the  ML  estimates  for  both  the  Weibull  and  K-  distribution  involve  iterative 
equations  manipulating  data  samples. 

Since  the  ML  estimates  are  readily  available  for  the  lognormal  distribution,  our 
discussion  focuses  on  the  Weibull  and  K-  distributions. 

A  few  estimation  methods  have  been  suggested  and  investigated.  Blacknell  (1994)  has 
compared  three  estimation  schemes  based  on  the  moment  method  for  K-distributions 
proposed  by  other  researchers.  Oliver  (1993)  has  used  the  Mellin  transform  to  form  an 
iterative  solution  to  achieve  the  approximate  ML  estimates.  Menon  (1963)  has  also 
suggested  a  way  to  obtain  estimates  for  the  Weibull  distribution.  In  this  report  we 
propose  two  estimation  schemes  named  as  no-bias  (NB)  I  and  II  estimates.  The  NB-I 
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estimation  scheme  estimates  parameters  with  no  bias  to  the  mean  and  to  the  variance 
in  the  log  domain  (£(z)  and  var(z) ),  whereas  the  NB-II  estimation  scheme  estimates 
parameters  with  no  bias  to  the  means  in  the  linear  domain  as  well  as  in  the  log  domain 
( E{x )  and  £(z)). 

4.3  NB-I  Scheme 

In  general,  for  a  perfect  two-parameter  distribution,  its  parameters  totally  determine 
the  mean  and  the  variance  and  vice  versa.  Intuitively,  as  we  are  interested  in  the  mean 
and  the  variance  of  a  distribution,  estimates  providing  no-bias  estimation  for  the  mean 
and  the  variance,  would  therefore  be  at  least  a  close  estimation,  if  not  the  optimum. 
Since  the  logarithmic  transform  is  non-linear,  and  often  we  are  interested  in  the 
distribution  presented  in  the  log  domain,  the  NB-I  scheme  simply  uses  the  first 
moment  e(z)  and  the  second  moment  e[z 2 )  to  estimate  parameters.  Therefore,  the 
scheme  is  no-bias  estimation  to  the  mean  and  variance  in  the  log  domain.  In  fact,  we 
have  seen  that  the  NB-I  estimates  are  also  the  ML  estimates  for  the  lognormal 
distribution. 

4.3.1  NB-I  Estimates  of  Weibull  Distribution 

Using  expressions  of  £(z)  and  var(z)  shown  in  Table  1,  the  NB-I  estimates  for  the 
Weibull  distribution  are, 


„  f  71 

c  =  exp  -y — j= 
V  -v  6 

*  71  1 

L  1. 

E{z)  "i 

sfc/(z)  J 

°  V6  Std(z) 

(26) 


It  can  be  proven  that  the  result  of  (26)  is  the  same  as  the  result  given  by  Menon  (1963) 
though  the  derivation  is  different.  The  physical  meaning  of  (26)  is  now  clarified,  that  is, 
the  mean  and  the  variance  in  the  dB  domain  of  the  Weibull  distribution  with 
parameters  b  and  c  are  equal  to  the  data  s  mean  E(z)  and  variance  var(z) , 
respectively. 

There  are  no  NB-I  estimates  for  the  K-distribution. 

4.4  NB-II  Scheme 

In  statistics  when  the  moment  method  is  used  for  parameter  estimation,  in  general  the 
higher  the  order  of  the  moment,  the  bigger  the  error  in  the  estimates.  The  NB-I  scheme 
uses  the  first  and  the  second  moments  to  derive  the  estimates.  Can  we  use  the  first 
moment  only  to  determine  the  estimates?  The  answer  is  yes.  In  order  to  estimate  two 
parameters,  we  may  use  E(x)  and  E{z) .  Since  the  transform  z  =  lOlogx  is  non-linear. 
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E(x )  and  E(z)  are  linearly  independent  though  both  of  them  are  the  first  moment  of 
the  data.  In  fact  E(x)  is  the  arithmetic  mean  and  E(z)  the  geometric  mean  (see 
Appendix  A  for  details). 

4.4.1  NB-II  estimates  of  Weibull  Distribution 

Using  the  expressions  of  E(x)  and  E(z)  shown  in  Table  1,  the  NB-II  estimates  for  the 
Weibull  distribution  are. 


In  r(l  +  a) +■ ^  +  ay  -  In  E(x)  =  0 


c  =  exp 


V 


E(z) 

k0a 


-Y 


(27) 


In  order  to  numerically  find  a  which  is  embedded  in  a  non-linear  equation,  iteration 
might  be  required.  But  this  iteration  is  different  from  the  iteration  involved  in  the  ML 
method.  This  iteration  does  not  need  the  involvement  of  the  sample  data  once  E(x)  and 
E{z)  are  ready.  Therefore  the  scheme  would  be  much  more  faster  than  the  scheme  in 
which  the  sample  data  are  involved  in  the  iteration. 

4.4.2  NB-II  Estimates  of  K-  Distribution 

For  the  K-  distribution,  obviously  the  best  estimation  of  x  is  E{x) ,  and  the  estimation 
of  the  a  can  be  found  from  the  expression  for  E{z)  shown  in  Table  2.  The  NB-II 
estimates  for  the  K-  distribution  are. 


x  =  e(x) 

'  if/^{a)-\aa  =  y- lnx  h — —  e(z) 

Similarly  parameter  a  can  be  numerically  determined  by  solving  the  non-linear 
equation  shown  in  (28),  which  is  again  much  faster  than  solving  the  non-linear 
equation  shown  in  (25)  for  the  ML  estimates  where  the  manipulation  of  the  sample 
data  is  involved  in  the  iteration. 

The  NB-II  scheme  for  the  K-  distribution  is  not  new,  and  a  similar  scheme  has  been 
studied  and  reported  by  Raghavan  (1991)  and  later  examined  by  Blacknell  (1994), 
though  the  derivation  is  different.  After  compared  three  non-ML  estimate  schemes, 
Blacknell  (1994)  concluded  that  the  NB-II  (though  Blacknell  did  not  use  such  a  name 
for  the  scheme  as  its  derivation  was  different  then)  was  the  best. 
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5.  Comparisons 


In  this  section  we  compare  different  estimation  schemes  using  both  the  RAR  and  SAR 
data. 

5.1  Low  Resolution  RAR  Data  of  Homogenous  Clutter 


Data  collected  by  the  multi-channel  airborne  radar  measurement  (MCARM)  system,  in 
particular  the  data  of  the  sum  analogue  channel,  extracted  from  Flight  5,  Circle  d. 
Acquisition  575,  are  used  for  this  case  study.  The  data  are  of  L-band  VV  low  resolution 
(the  ground  patch  size  is  about  8xl05m2)  and  represent  the  clutter  of  statistically 
homogenous  farmland.  The  total  sample  number  is  762,  which  seems  to  be  small  for 
the  distribution  analysis,  but  this  is  the  maximum  number  we  can  extract  from  the  data 
set.  The  data  need  pre-processing,  and  the  details  are  given  in  Appendix  B.  Since  the 
processing  gains  in  various  stages  are  unknown,  the  mean  of  the  samples  is  normalised 
to  1  (OdB). 

The  estimates  calculated  from  the  ML,  NB-I  and  NB-II  estimation  schemes  for  the 
Weibull  and  K-  distributions,  respectively,  are  given  in  Table  5.  To  evaluate  the 
accuracy  of  the  NB-I  and  NB-II  estimation  schemes,  we  can  use  the  ML  estimates  as  the 
benchmark,  as  the  ML  estimates  are  optimal.  According  to  the  estimated  parameters, 
the  distribution  is  very  close  to  the  Rayleigh  distribution,  which  means  the 
corresponding  parameter  a  of  the  K-distribution  should  approach  infinity.  The  values 
of  a  calculated  by  the  ML  and  NB-II  schemes,  respectively,  seem  to  be  significantly 
different,  but  it  is  merely  the  calculations  converging  to  different  points.  For  numerical 
calculation,  all  a >25  can  be  regarded  as  the  same  as  this  will  become  clear 

when  we  look  at  the  distribution  curves. 

Table  5:  ML  NB-I  and  NB-II  estimates  of  the  MCARM  RAR  (L  band  VV  polarisation)  clutter 
data  of  farmland. 


Distribution 

ML  estimates 

No-bias  estimates 

Scheme  I 

Scheme  II 

Weibull 

a  =  0.9916 

c  =  0.9963 

a  =  0.9359 

c  =  1.0118 

a  =  0.9738 

c  =  0.9889 

K- 

3c  =  1.000 

6  =  25.5876 

N/A 

3c  =  1.0000 

a  =  72.0000 

The  pdfs  and  cdfs  of  the  Weibull  distribution  using  the  ML,  NB-I  and  NB-II  estimates 
are  plotted  in  Figure  6,  while  the  pdfs  and  cdfs  of  the  K-distribution  using  the  ML  and 
NB-II  estimates  are  plotted  in  Figure  7.  In  order  to  examine  the  details  of  the 
distribution  at  the  ends,  the  pdfs  on  the  log  scale  are  also  plotted.  The  cdfs  are  also 
plotted  on  the  log  scale,  since  we  are  more  interested  in  the  high  end. 
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It  can  be  seen  that  the  difference  between  the  ML  estimates  and  the  NB-II  estimates  for 
both  the  Weibull  and  the  K-  distributions,  respectively,  is  insignificant,  and  the 
corresponding  pdfs  and  cdfs  are  nearly  identical.  The  NB-I  estimates  are  also  good,  but 
not  as  good  as  the  NB-II  estimates. 

The  comparison  between  the  Weibull  and  the  K-  distributions  are  shown  in  Figure  8. 
The  K-  distribution  converges  a  little  bit  slower.  If  we  view  Figure  6  and  Figure  7 
carefully,  the  K-distribution  seems  a  little  better  to  represent  the  data  distribution 
especially  at  the  higher  end. 
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Figure  6:  Comparison  of  parameter  estimation  schemes  for  the  Weibull  distribution  using  the 
MCARM  RAR  L-band  VV  clutter  data  of  homogenous  farmland. 
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5.2  AirSAR  Single-Look  L-band  VV  Data  of  Homogenous  Clutter 

SAR  synthesises  the  length  of  radar  antenna  to  improve  the  azimuth  resolution.  If  a 
RAR  system  had  the  same  frequency,  polarisation,  resolution  and  SNR  as  that  of  a  SAR 
system,  the  collected  data  would  be  statistically  the  same.  In  this  case  study  we  analyse 
a  block  of  the  AirSAR  single-look  L-band  W  data,  representing  the  clutter  of  a 
statistically  homogenous  open  woodland  of  native  eucalyptus.  The  incidence  angle  of 
the  block  is  about  35-45°,  well  within  the  plateau  region  of  the  grazing  angle  for  forests. 
The  clutter  can  be  considered  as  independent  of  incidence  angle  (Dong,  2003).  The 
variation  is  therefore  purely  spatial.  The  block  consists  of  375x336  pixels,  enough 
sample  numbers  for  deriving  reliable  statistics.  The  pixel  size  is  about  3.6x6.7  m2.  The 
data  was  processed  by  JPL  back  in  1993  as  non-standard  single-look  data  (contrast  to 
the  standard  multi-look  data). 

The  parameters  estimated  by  the  ML  and  NB-I  and  NB-II  schemes  are  given  in  Table  6. 
The  corresponding  cdfs  and  pdfs  for  the  Weibull  and  K-  distributions  are  shown  in 
Figure  9  and  Figure  10,  respectively.  The  distribution  is  again  very  close  to  the 
Rayleigh  distribution.  The  differences  among  all  the  estimates  for  the  Weibull  and  K- 
distributions,  respectively,  are  very  small  and  the  pdfs  and  cdfs  plotted  using  these 
estimates  are  nearly  identical. 

The  comparison  between  the  Weibull  and  the  K-  distributions  are  shown  in  Figure  11. 
As  can  be  expected,  the  K-  distribution  converges  a  little  bit  slower. 


Table  6:  Comparison  of  parameter  estimation  schemes  for  AirSAR  single-look  L-VV  data. 


Distribution 

ML  estimates 

No-bias  estimates 

Scheme  I 

Scheme  II 

Weibull 

a  =  1.0243 

c  =  2.4536 

a  =  1.0097 

c  =  2.4947 

a  =  1.0197 

c  =  2.4583 

K- 

x  =  0.4030 

&  =  23.7554 

N/A 

x  =  0.4030 

a  =  25.3648 
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5.3  AirSAR  Single-Look  C-band  VV  Data  of  Heterogeneous  Clutter 

In  this  case  study  we  analyse  a  block  of  the  AirSAR  single-look  C-band  VV  data. 
Instead  of  selecting  homogenous  clutter,  we  deliberately  select  a  block  containing  open 
eucalypt  woodland,  farmland  and  scattered  farmhouses.  The  block  consists  of  300x400 
pixels,  enough  sample  numbers  for  deriving  reliable  statistics.  The  size  of  the  pixel  is 
the  same  as  die  AirSAR  L-band  data. 

The  parameters  calculated  using  the  ML,  NB-I  and  NB-II  schemes  are  given  in  Table  7. 
Since  the  clutter  is  not  homogenous,  the  distribution  becomes  broader.  This  might  be 
considered  as  the  sum  of  the  two  Rayleigh  distributions  with  a  few  dB  mean  difference. 
The  corresponding  pdfs  and  cdfs  for  the  Weibull  and  K-  distributions,  respectively,  are 
shown  in  Figure  12  and  Figure  13.  For  the  Weibull  distribution  the  NB-I  scheme  does 
not  produce  as  good  estimates  as  the  NB-II  scheme  in  this  case.  The  difference  between 
the  cdfs  of  the  ML  and  NB-I  estimates  grows  to  about  ldB  when  the  cdf  reaches  0.9999. 
The  difference  between  the  ML  and  NB-II  estimates  for  both  the  Weibull  and  the  K- 
distributions,  respectively,  however  is  again  insignificant,  and  the  corresponding  pdfs 
and  cdfs  are  nearly  identical. 

The  comparison  between  the  Weibull  and  the  K-  distributions  are  shown  in  Figure  14. 
Since  the  distribution  widens,  but  the  parameter  a  of  the  Weibull  distribution  is  still  in 
the  range  of  1  <  a  <  2 ,  the  convergence  of  the  K-distribution  is  slower  than  that  of  the 
Weibull  distribution. 


Table  7:  Comparison  of  parameter  estimation  schemes  for  the  AirSAR  single-look  C-VV  data  of 
heterogeneous  clutter. 


Distribution 

ML  estimates 

No-bias  estimates 

Scheme  I 

Scheme  II 

Weibull 

d  =  1.2009 

c  =  2.8894 

a  =  1.1312 

c  =  3.1324 

a  =  1.1855 

c  =  2.8953 

K- 

x  =  0.3100 

a  =  2.7389 

N/A 

3c  =  0.3100 

a  =  2.7035 
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5.4  High  Resolution  Radar  X-band  VV  Data 

The  data  was  collected  by  a  surface-mounted  X-band  (centre  frequency  of  9750  MHz) 
W  polarised  radar  with  a  range  resolution  of  0.3  m  and  PRF  of  1000  Hz  at  a  grazing 
angle  of  approximate  zero  degrees  looking  at  the  sea  surface.  The  original  data  set 
consists  of  6000  pulses  by  830  range  bins.  The  range  of  the  first  bin  is  5425.5  m.  The 
change  in  range  for  other  range  bins  is  therefore  insignificant,  so  the  received  power 
can  be  considered  as  independent  of  range,  and  its  variation  can  be  considered  to  be 
purely  due  to  the  spatial  and  temporal  distribution.  The  data  collected  in  the 
corresponding  range  bins  for  the  consecutive  pulses  are  often  highly  correlated  and 
cannot  be  considered  as  independent  samples.  The  correlation  between  pulse  i  and 
pulse  j  is  defined  as  (Proakis  and  Manolakis,  1996), 


Pij  = 


Q 

X«i*« 

k= 1 


jk 


9  *  9  * 

X Y.ujkujk 

k=l  k= 1 


1/2 


(29) 


where  *  denotes  complex  conjugate,  uik  is  the  radar's  ith  pulse  measurement  in  range 
bin  k ,  and  q  is  the  total  number  of  range  bins. 

Figure  15  shows  the  correlation  between  the  data  collected  in  pulse  1  and  the  following 
pulses.  It  is  safe  to  say  that  when  the  interval  of  pulses  is  greater  than  about  300  (300 
milliseconds),  the  data  become  effectively  uncorrelated.  In  our  calculation,  an  interval 
of  300  pulses  was  used,  so  the  data  collected  for  pulse  1,  301,  601,  ...,  and  5701  were 
treated  as  independent  samples  given  a  total  sample  number  to  be 
6000/300x830  =  16600 .  It  is  worth  noting  that  the  variation  of  the  selected  data  is  both 
spatial  and  temporal. 

The  mean  of  the  data  was  normalised  to  1  (OdB).  The  distribution  was  found 
significantly  different  from  the  previous  cases  and  close  to  a  lognormal  distribution. 
This  might  be  due  to  the  illumination  geometry  of  high  resolution  and  low  grazing 
angle.  A  high  resolution  means  that  an  illuminated  patch  does  not  contain  many 
distributed  but  a  few  discrete  scatterers  while  a  low  grazing  angle  creates  effects  of 
shadowing  and  multipath  propagation.  Therefore  the  criterion  of  Rayleigh  distribution 
does  not  hold.  The  pdfs  of  the  data  as  well  as  the  ML  lognormal,  Weibull  and  K- 
distribution  curve  fitting  are  shown  in  Figure  16.  It  can  be  seen  that  the  pdf  of  the  data 
is  not  quite  a  lognormal  distribution,  but  among  the  three  distributions,  the  lognormal 
distribution  is  the  best  to  fit  the  data.  The  ML  lognormal  distribution  is  even  a  little 
conservative  at  the  high  end  (the  disagreement  in  the  low  end  might  be  due  to  the 
noise  level  of  the  radar)  which  can  be  further  identified  by  viewing  the  cdfs  of  the  data 
and  the  ML  lognormal  curve  fitting  as  shown  in  Figure  17. 
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Figure  17:  Cdfs  of  the  data  and  the  ML  estimated  lognormal  curve  fitting.  Since  the  true  data 
distribution  is  not  quite  a  lognormal  distribution,  the  ML  estimates  are  a  little  conservative  at 
the  high  end. 


The  ML  estimated  parameters  for  the  three  distribution  functions  are  given  in  Table  8. 
However,  it  has  to  be  pointed  out  that  pdf  parameters  of  sea  clutter  are  highly 
dependent  on  radar  frequency,  polarisation,  grazing  angle,  as  well  as  weather 
conditions. 

Table  8:  The  ML  estimated  parameters  for  the  lognormal,  Weibull  and  K-  distributions  for  the 
approximate  zero  degree  grazing  angle,  high  resolution  X-band  VV  sea  clutter  data  (the  mean  of 
tire  data  has  been  normalised  to  1). 


Distribution  function 

Parameters 

Lognormal 

//  =  -1.4892;  5  =  1.7207 

Weibull 

b  =  0.5799;  c  =  1.4453 

K- 

x  =  1.0;  a  =  0.8080 

6.  Conclusions 


The  Weibull,  K-  and  lognormal  distributions,  widely  used  to  represent  radar  clutter 
distributions,  have  been  studied. 

Since  clutter  is  often  expressed  in  the  dB,  its  pdf  and  cdf  are  also  plotted  on  the  dB 
scale,  the  statistical  properties  including  mean  and  variance  for  the  Weibull  and  K- 
distributions  in  the  log  domain  have  been  derived.  These  values  are  later  used  for 
parameter  estimations. 


If  the  clutter  is  modelled  as  a  Weibull  distribution,  the  shape  parameter  a  is  generally 
in  the  range  of  l<a<2.  Corresponding  to  this  range,  the  lognormal  distribution 
converges  slowest  among  the  three  distributions,  followed  by  the  K-distribution,  and 
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the  Weibull  distribution  converges  fastest.  This  finding  contradicts  a  general  claim  of 
Billingsley  (2002)  that  the  K-distribution  converges  faster  than  the  Weibull  distribution. 
In  fact  his  claim  is  true  only  when  the  distributions  correspond  to  a  >  2 . 

However,  the  value  of  a  generally  becomes  larger  if  the  grazing  angle  is  low  and  the 
resolution  is  high.  The  value  of  a  ranging  from  1.8  to  5.4  for  data^of  various  terram 
types  including  urban  and  forests,  with  the  grazing  angle  less  than  2  and  resolution  m 
the  order  of  10V,  has  been  reported  (Billingsley,  2002).  If  the  resolution  reduces  to 
106m2,  the  value  of  a  becomes  1.1  to  2.8  (Billingsley,  2002). 

According  to  both  the  RAR  and  SAR  data  analysed,  visually  the  K-distribution  seems 
marginally  better  than  the  Weibull  distribution  in  describing  the  clutter  distribution 
especially  at  the  high  end  of  the  distribution.  We  did  not  use  the  standard  or  modified 
statistical  tests  such  as  the  chi-square  or  modified  chi-square  tests  (Antipov,  1998, 
Chan,  1990)  to  numerically  test  which  distribution  is  better.  Because  we  are  only 
concerned  about  the  distribution  at  the  high  end,  the  error  occurring  in  other  parts  of 
the  distribution  should  not  be  significant  for  this  programme. 

Theoretically,  if  the  surface  is  homogenous  and  contains  many  scatterers,  the 
distribution  of  the  radar  clutter  is  Rayleigh.  Both  the  RAR  and  SAR  data  analysis 
support  this  theory.  The  radar  resolution  of  the  RAR  data  is  of  the  order  of  tens  of 
thousand  square  metres  whereas  the  resolution  of  the  AirSAR  single-look  data  is  of  tire 
order  of  tens  of  square  metres,  but  both  data  approach  the  Rayleigh  distribution  for 
homogenous  clutter.  For  heterogenous  areas,  the  clutter  distribution  is  broader  than 
the  Rayleigh  distribution.  This  can  be  considered  as  the  sum  of  two  or  more  Rayleigh 
distributions  with  slightly  different  means.  The  distribution  of  sea  clutter  of  very  fine 
resolution  and  low  grazing  angle  is  approximately  lognormal,  accordmg  to  the  data 
analysed.  At  low  grazing  angle,  the  shadowing  and  multipath  propagation  effect  may 
dominate  so  the  Rayleigh  criterion  does  not  hold. 

Parameter  estimation  is  an  important  issue  and  the  ML  estimates  are  the  optimal. 
Previous  research  concentrated  on  using  the  Lagrange  method  to  find  the  ML 
estimates.  For  the  K-distribution  because  the  derivative  of  the  embedded  Bessel 
function  cannot  be  expressed  in  a  closed  form,  it  is  very  difficult  to  find  the  ML 
estimates  using  the  Lagrange  method,  if  not  impossible.  This  report  has  proposed  a 
direct  method  to  find  the  ML  estimates  for  the  K-distribution.  The  direct  method  does 
not  require  the  derivative  of  the  Bessel  function. 

The  ML  method  provides  the  optimal  parameter  estimates.  Unfortunately,  except  for 
the  lognormal  distribution,  the  ML  estimates  for  both  the  Weibull  and  K-  distribution 
involve  iterative  computation  manipulating  sample  data.  Faster  estimations  are 
desirable  in  real-time  data  processing.  This  report  has  proposed  two  no-bias  estimation 
schemes,  namely  NB-I  estimates  and  NB-II  estimates.  The  NB-I  scheme  uses  the  mean 
and  the  variance  of  data  in  the  log  domain  to  calculate  the  parameters  (so  the  estimated 
mean  and  variance  equal  the  mean  and  variance  of  the  data  in  the  log  domam),  while 
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the  NB-II  scheme  uses  the  mean  of  the  data  in  the  linear  domain  (arithmetic  mean)  and 
the  mean  of  the  data  in  the  log  domain  (geometric  mean)  to  calculate  the  parameters 
(so  the  estimated  means  equal  the  means  of  die  data  in  the  linear  and  log  domains, 
respectively).  It  has  been  found  that  NB-II  performs  better,  as  the  estimation  only  uses 
the  first  moment  of  the  data.  The  differences  between  the  NB-II  estimates  and  die  ML 
estimates  are  insignificant,  and  the  corresponding  pdfs  and  cdfs  are  nearly  identical  for 
all  cases  studied.  The  accuracy  of  the  NB-I  is  not  as  good  as  that  of  NB-II,  though  the 
difference  is  still  small. 
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Appendix  A:  Statistics  of  Weibull  and  K- 
Distributions  in  Log  Domain 

In  order  to  maintain  the  physical  meaning,  generally  radar  signals  such  as  power,  RCS 
and  reflectivity  are  manipulated  in  the  linear  domain,  though  their  values  are  often 
expressed  in  dB.  For  instance,  the  average  power  of  0  dBw  (1  w)  and  10  dBw  (10  w)  is 
7.4  dBw  (5.5  w)  but  not  5  dBw  (3.16  w),  unless  specified.  However,  if  we  define  the 
manipulation  domain,  we  may  accept  7.4  dBw  as  the  arithmetic  mean  and  5dBw  as  the 
geometric  mean.  Without  loss  of  generality,  we  define  the  mean  or  more  precisely  the 
arithmetic  mean  as. 


1  « 

xarith  ~  ^lxi 
«/= 1 


(Al) 


Similarly  we  may  define  the  geometric  mean  as. 


( »  Y/n 

xgecm  =1 fl*/J  */>0  (A2) 

When  the  value  of  xt  is  expressed  in  dB,  (A2)  is  equivalent  to  as, 

z=-±z,  (A3) 

«/= 1 

where  z  =  10 log*  . 

Therefore  we  may  consider  the  manipulation  in  the  linear  domain  as  the  arithmetic 
manipulation  and  the  manipulation  in  the  log  domain  as  the  geometric  manipulation. 
Since  radar  signals  are  often  expressed  in  the  dB  domain,  the  geometric  manipulation 
sometimes  may  have  advantages.  This  Appendix  derives  statistics  for  the  Weibull  and 
the  K-  distributions  in  the  log  domain.  The  statistics  of  the  lognormal  distribution  in 
the  log  domain  is  simply  the  statistics  of  a  normal  distribution. 

A.l.  Statistics  of  Weibull  Distribution  in  Log  Domain 

The  pdf  of  the  Weibull  distribution  in  the  natural  logarithmic  domain  is, 

p  0)  =  bc{ey  J  expf- c(ey  f  1  -oo  <  y  <  oo  ( A4) 


where  y  =  lnx  . 
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To  find  the  mean  and  variance  of  the  distribution  requires  the  evaluation  of  the  first 
and  the  second  moment  of  the  distribution. 


E(y)=  \yp(y)<fy 


E(y2)= 

-co 

Let  t  =  c(ey  J ,  we  have 


(A5) 

(A6) 


(AT) 


WJ  w 

In 2  c  +  J  In 2 1  e-'  dt  -  2  In  c  J  In  t  e~l  dt 


(A8) 


According  to  Gradshteyn  and  Ryzhik  (1980),  and  Mathemetica,  Version  4  (Wolfram 
Research,  2000)  we  can  evaluate  the  special  integrals  in  (A 7)  and  (A8)  as. 


| In  te  1  dt  =  -v/(0)  (l)  =  Y  (A9) 


where  y/^ix)  is  the  Polygamma  function,  defined  as  the  nth  derivative  of  the  logarithm 

jr/l+1 

of  the  Gamma  function,  — — (lnr(x))  and  y  is  Euler's  Gamma  constant. 

dx 


fin 2  te  ’dt  =  y 1  +—■ 
o  6 

Finally  we  have, 

E(y)  =  ~(y  +  lnc) 
b 

var(y  )  =  1 )- E2  (y)  = 

6  b 

In  the  dB  domain,  the  above  expressions  become. 


(A10) 


(All) 

(A12) 


38 


DSTO-RR-0274 


e(z)  =  -^-(y  +  Inc) 


var(z)  = 


AoV 

6b2 


(A13) 

(A14) 


where  k0  =  10/lnlO 

A.2.  Statistics  of  K-  Distribution  in  Log  Domain 

The  pdf  of  the  K-  distribution  in  the  natural  logarithmic  domain  is. 


(A15) 


Letting  t  =  2(01  ey  /  x )'  ,  we  have 


E^’r{a]2.-d(2'a,+'n^KU,)dl 

(A16) 

E^-r(al°4(2b,+bto)2,°K‘-'{^ 

(A17) 

Evaluations  of  the  special  integrals  in  (A16)  and  (A17)  include  (Gradshteyn  and 
Ryzhik,  1980,  Mathemetica,  Version  4, 1998), 

<jtaKa_1(t)dt  =  2“-1r(or) 

(A18) 

0 


]ta\nt  KaA  (t)dt  =  2a~2  r(or)[- /  +  21n2  +  /0)  (or)]  ( A19) 

0 


]ta  In2  tKa^(t)dt 
0 

=  2a_4r(a^2r2  +—^ — 8rln2  +  81n22-4(r-21n2V(0)(a)+2(^(0)(a))2  +2^(1)(ar) 

(A20) 


Finally  we  have. 
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£(y)=  -y +y/^(pc)-  Ina  H-lnx 

var(y)  =  ^(1)(«)+^~ 

o 

In  the  dB  domain,  the  above  expressions  become  as, 
E(z)  =  k0  (-  y  +  ^(0)  (°0 -  In  a  +  In  * ) 

var(z)=  k02  iy{l){a)+~- 


(A21) 

(A22) 

(A23) 

(A24) 


We  have  proven  in  Section  2.2  that  the  K-  distribution  with  or  —  1/2  is  identical  to  the 
Weibull  distribution  with  6  =  1/2  and  c  =  ( Hx )1/2 .  In  fact, 

y/(0)(l/2)  =  -r-21n2  (A25) 

^(1)(l/2)  =  7T2  /2  (A26) 

We  can  see  that  the  statistics  of  the  Weibull  and  K-  distributions  are  indeed  the  same 
with  these  special  conditions. 
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Appendix  B:  RAR  Data  Pre-processing  for  Spatial 

Distribution  Analysis 


According  to  the  radar  equation,  the  received  power  from  terrain  is, 

PavgGtGrX2<TQFUg 
r  {An  )3i?4I 


(Bl) 


where  Pavg  -  average  transmitted  power,  G,  and  Gr  -  gains  of  radar  transmitting  and 
receiving  antenna  in  the  direction,  A  -  wavelength,  Ag  -  resolvable  area  of  the  ground, 
R  -  range  of  the  ground  patch,  <r0  -  reflectivity,  F*  -  the  pattern  propagation  factor 
and  L  the  total  loss  including  propagation  loss  and  system  loss. 

In  order  to  analyse  the  spatial  distribution  of  clutter,  cr0F4,  we  need  to  pre-process  the 
received  power  in  such  a  way  to  extract  cr0F* .  This  is  done  in  two  steps.  First  we 
calculate  the  clutter  that  is  only  a  function  of  grazing  angle.  Secondly  we  truncate  those 
data  which  are  highly  correlated  to  the  grazing  angle.  The  distribution  of  the  remaining 
data,  which  are  independent  of  (almost  uncorrelated  to)  the  grazing  angle,  therefore 
can  be  considered  as  the  spatial  distribution. 

For  a  given  range,  the  area  illuminated  by  the  3dB  main  beamwidth  can  be  written  as, 
Ag  =  rAr<f>3dB  (B2) 


where  r  »  R  cos  9 ,  9  is  the  depression  angle;  A r  = 


rc 


2  cos  P 


,  P  is  the  grazing  angle  and  r 


is  the  time  interval  between  two  consecutive  range  bins.  For  the  MCARM  system,  the 
height  of  the  aircraft  is  about  3500  m,  r  =  0.8 jus .  With  these  parameters,  and  noting  the 
grazing  angle  is  very  low,  so  cos 9  ~  cos/?,  the  ground  resolution  is  in  the  order  of 
8  x  105  m1 ,  a  very  coarse  resolution  for  the  analogue  sum  channel. 

Assuming  the  clutter  in  the  3dB  main  beamwidth  area  to  be  constant,  we  have. 


< TnF 


PrSdB  °c-VL  Gt(9J)Gr(9JW 


R 


(B3) 


where  Pr3dB  is  the  power  received  from  the  area  illuminated  by  the  main  beam  of  3dB 
beamwidth.  When  radar  looks  at  boardside,  the  Doppler  resolution  is  in  the  azimuth 
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direction.  Therefore,  after  the  received  power  being  transformed  into  the  Doppler 
frequency  domain,  the  corresponding  component  of  Pr3dB  can  be  determined.  The 
Doppler  frequency  shift  due  to  the  movement  of  vegetation  at  L-band  is  relatively 
small  (Billingsley,  2002),  and  the  main  contribution  to  the  Doppler  frequency  shift  is 
therefore  due  to  the  geometry. 

!•«  -  £w  (M> 

/=/, 

where  P(f)  is  the  power  component  of  the  frequency  /  and 


2v 

fx  =  —  cos  6  cos] 


n  ,  fa  dB  c 

- i - Oq 

2  2 


(B5) 


2v  i 

f0  = — cos  9  cos 
J  2  A  1 


(B6) 


where  v  is  the  velocity  of  the  platform,  9  the  depression  angle,  S0  the  drift  angle  of 
the  platform  and  fadB  the  azimuth  angle  of  3dB  beamwidth. 


We  then  have  the  final  expression. 


a0FA  x 


PrldB^ 

G2{6) 


(B  7) 


where 

g 2  (9)  =  Gt  ( 9 ,  <f>)Gr  ( 9 ,  fad<f>  (B8) 


The  above  integration  is  along  the  contour  of  constant-elevation,  which  also  assumes 
the  range  on  the  contour  of  constant-elevation  is  the  same.  Rigorously  speaking,  the 
contour  of  constant  elevation  angle  of  a  tilt  antenna  array  on  the  ground  is  nearly  a 
straight  line  whereas  the  contour  of  the  equal-range  is  a  half-circle.  However,  we  can 
consider  the  two  contours  to  be  coincident  for  a  small  angle  of  faB . 

The  result  of  (B 7),  the  clutter  as  a  function  of  grazing  angle  is  shown  in  Figure  Bl. 
Apart  from  the  notch  which  corresponds  to  the  clutter  of  bay  water,  the  clutter  shown 
in  the  figure  is  of  low  relief  farmland.  The  clutter  of  farmland  seems  to  be  independent 
of  grazing  angle  when  the  angle  is  grater  than  5°,  and  the  variation  can  be  considered 
as  the  random  process  of  clutter  in  space.  As  shown  in  Figure  B2  range  bins  200-316 
(because  of  the  transmitter  signal  leakage  in  range  bin  68,  signals  in  range  bins  of  1-200 
are  not  used)  satisfy  this  criterion  of  9  >  5° .  This  provides  a  total  of  127  samples,  too 
less  for  the  spatial  distribution  analysis. 
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- Medium  PRF 


Figure  Bl:  Clutter  as  a  function  of  grazing  angle.  The  notch  corresponds  to  clutter  of  bay  water 
and  the  others  to  farmland. 


Figure  B2:  Grazing  angle  and  depression  angle  as  functions  of  range  bins. 

The  pulse  number  in  the  CPI  process  of  the  medium  PRF  MCARM  data  is  128. 
Unfortunately  the  corresponding  samples  collected  in  different  pulses  are  generally 
highly  correlated,  and  therefore  are  not  independent  samples.  Figure  B3  shows  the 
correlation  between  the  corresponding  samples  collected  from  pulse  1  and  pulse  k 
( k  =  1,2,  •••,128 ).  Correlations  among  other  CPI  pulses  have  the  similar  pattern.  It  can  be 
seen  that  in  general  the  samples  become  independent  if  the  interval  between  the  CPI 
pulses  are  more  than  20.  Therefore  we  have  selected  samples  collected  from  pulses  1, 
25,  50,  75,  100  and  125  for  the  distribution  analysis.  The  total  sample  number  is 
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therefore  127x6  =  762 .  The  number  is  still  not  sufficient  for  distribution  analysis,  but 
this  is  the  maximum  number  we  can  extract  from  one  medium  PRF  MCARM  data  set. 


Figure  B3:  Signal  correlation  between  the  1st  CPI  pulse  and  the  kth  CPI  pulse  (k  -  12, •••,128). 
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